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Abstract: The parametric resonance responsible for preheating after inflation will end 
when self-interactions of the resonating field and interactions of this field with secondary 
degrees of freedom become important. In many cases, the effect may be quantified in 
terms of an effective mass and the resulting shifting out of the spectrum of the strongest 
resonance band. In certain curvaton models, such thermal blocking can even occur before 
preheating has begun, delaying or even preventing the decay of the curvaton. We investigate 
numerically to what extent this thermal blocking is realised in a specific scenario, and 
whether the effective mass is well approximated by the perturbative leading order thermal 
mass. We find that the qualitative behaviour is well reproduced in this approximation, and 
that the end of preheating can be confidently estimated. 
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1 Introduction 

At the end of inflation, the Universe was reheated through the transfer of energy from 
the inflaton field to a number of coupled scalar, vector and fermion helds. This occurred 
through a combination of perturbative decay and non-perturbative preheating. Primary 
particles created in such processes subsequently decayed into other secondary degrees of 
freedom, and the whole system ultimately thermalised to produce a radiation-dominated 
Universe. 

Resonant preheating occurs because the post-inflationary oscillations of the inflaton are 
in resonance with certain momentum modes of the primary field(s) [1]. The final particle 
number depends on the strength and duration of the resonance. In general, there is an 
infinite tower of resonance bands, but the strength of the resonance drops significantly 
beyond the first band. In the simplest case, a first-band resonance occurs if there is a 
momentum mode k satisfying 


k^ + AfJj — mfnf, (1.1) 

where ruinf is the frequency of inflaton oscillation. Meg is the effective mass of the held 
coupled to the inhaton; it changes in time due to both interactions with other helds and 
self-interactions. This paper deals with the blocking of the resonance due to the back¬ 
ground temperature causing an effective thermal mass to be generated (termed ’’thermal 
blocking”). This could happen either as the system heats up, such that a resonance that is 
effective at the beginning of the preheating process shuts off when the temperature (or out- 
of-equilibrium spectrum) is established, or due to an already-present thermal background, 
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such as in the curvaton model. We do not discuss in this paper the backreaction of the 
resonant field on the oscillating field, which modifies the oscillation frequency and can also 
end preheating. 

The same physical process can occur in curvaton models [2] , where a scalar field distinct 
from the inflaton is present, but has negligible energy density during inflation and decays 
long after the end of inflation. Before reheating, the curvaton comes to almost dominate 
energy density of the Universe, and upon reheating, its isocurvature perturbations become 
the adiabatic curvature perturbations observed in the CMB. Similarly to the inflaton case, 
the curvaton can decay either perturbatively or through non-perturbative preheating. In 
contrast to the inflaton case, there generically exists a thermal background at the start of 
curvaton decay, formed by the inflaton’s decay products. The process of curvaton reheating 
has been the subject of recent attention [3], particularly the report that preheating in the 
curvaton model into Higgs particles can be substantially delayed by the presence of the 
thermal background [4] . This result is an integral component of the minimal curvaton-higgs 
(MCH) model, in which the curvaton is coupled directly to the Standard Model Higgs boson 
[5]. Without this blocking, the model would become unviable, unless the Higgs-curvaton 
coupling were tuned to be tiny, because reheating would occur while the curvaton was too 
subdominant. Thus, the viability of the MCH model depends on whether the analytical 
calculations in [4] hold up beyond the simple LO perturbative estimate of the thermal 
mass. 

The onset and end of non-perturbative instabilities has been studied numerically for 
some time in the context of resonant preheating (see for instance [6]), tachyonic preheating 
(see [7]), self-resonance (see [8]), and also for heavy-ion collisions (see [9]). Because the 
particle numbers are often very large, a classical-statistical approach is usually adopted, 
but also quantum dynamics has been studied in this context [10, 11]. In most cases, the 
main emphasis has been on the kinetic thermalisation phase and the back-reaction on the 
field driving the instabilities, rather than in terms of a thermal mass. For a general review 
of nonperturbative reheating, see [12]. 

The two aims of this paper are to determine the conditions under which the thermal 
blocking effect is present; and to determine the extent to which thermal blocking is well 
described by an effective leading order thermal mass. 

Our paper is organised as follows. In section 2 we introduce our model and our analysis 
of preheating, including the concept of the thermal mass in the leading order (LO) approx¬ 
imation. We also discuss the timescale of the various physical processes. In Section 3 
we present our numerical results for thermal blocking, showing both the results in the 
LO approximation and once dynamical secondary fields are introduced. We conclude in 
section 4. 

2 Thermal blocking of preheating 

Our numerical investigation uses the curvaton model as an example. This is because 
the thermal background is produced well before curvaton preheating, by the decay of the 
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inflaton. Thus, it is much easier to model than a pure inflaton model, where the thermal 
background is produced during preheating. 


2.1 Model 

We consider a very simplified model of the Standard Model Higgs field (p coupled to both 
the curvaton scalar a and Nf additional light degrees of freedom i = These 

light fields represent the fermions and gauge fields of the Standard Model, in terms of decay 
channels and thermal back-reaction on the Higgs field. We take them to be massless. 

We will assume that the Standard Model degrees of freedom {cj) and the are initially 
in thermal equilibrium at a temperature T, and that a is homogeneous and oscillating. The 
field a is coupled only to the Standard Model Higgs. Due to this coupling, there is potential 
for resonant production of Higgs excitations. The additional fields act as secondary fields; 
they are not coupled to the oscillating field a. We neglect the expansion of the Universe. 

The Standard Model Higgs field is a 2-component complex scalar field given by 


1 / </>i + \ 

^/2 yh + j ' 


( 2 . 1 ) 


and thus the action is 


S = - d^x 




( 2 . 2 ) 


with the known constants 


m 


g = 


H 


126 ^ GeV^ 


A 


= v^ = 246^ GeV^ 


(2.3) 


and an unknown coupling g. The function = cr‘^{t) is a priori unknown; for simplicity 
we take it to be 


cj^(t) = CTg cos^ {niat ), (2.4) 

with some new unknown constants do, and rua- Note that, in a fully dynamical setup, the 
field a would also be dynamical and interact with its environment (we do not implement 
this here). 

The equation of motion for the Higgs field is then 




+ 2X(p^ (p + g‘^a‘^{t) + 

i 


(j){x,t) = 0, 


(2.5) 


and for the secondary fields. 


-|- 2h?4)^(1) 


^i{x,t) = 0. 


( 2 . 6 ) 


Given g, do, rua and h, these can be solved numerically. We have chosen the secondary 
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fields to be massless, which is approximately true for reheating temperatures beyond the 
electroweak scale (except for the top quark). Because the top Yukawa coupling is of order 
one, and the electroweak gauge coupling is not much smaller than unity, the effective 
Higgs-SM coupling h should also be taken to be of order one. 


2.2 Resonant preheating 

In our model, both primary and secondary fields start out in thermal equilib¬ 

rium. Thus, the Fourier components of each field and its conjugate momentum tt* = dt^ 
(in continuum notation), which are given by 

obey^ 

(4(k)(/)i(k)) = —, (7r|(k)7ri(k)) = nkWk, (2.8) 

Wk 

with -|- |k|^. Particle numbers are given by the Bose-Einstein distribution 

_ 1 ■ P-f) 

Similar expressions hold for the fields. 

The equation of motion for the Higgs field in momentum space is (including only local 
self-energy components; see section 2.3 for details) 




+ ikr- 


cos^' 


{m„t) (/)k = 0. 


( 2 . 10 ) 


A good choice for the initial state mass of the Higgs field is 


"^fnitial — + Af^(T) + Gq, 


( 2 . 11 ) 


with the assumption that this quantity is positive. As stated previously, are massless. 
The possibility of resonant preheating can best be seen by rewriting (2.10) into the Mathieu 
equation using 


Ak — 


|k|^ -g? + + M^{t) + gVp 

m'i 


/2 


2^2 


Q = 


g q'l 

4m^ 


T = m. 


,t, (2.12) 


to get 


+ (Ak - 2gcos(2T)) (/>k = 0. 


(2.13) 


^We do not include the “quantum half” (zero point energy). Doing so would require renormalisation of 
the mass, but otherwise results in qualitatively similar results. 
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This equation has resonant solutions around ^4^ = 1,4,9,corresponding to 


1^ ^ + - A((/»V) - M^(t) - gVg/2 

rri^ 

In the following, we take the fiducial values = 0.1, (jQlra„ = 1, and tuh/ rna = 0.3. Note, 
that because a is not a dynamical field, only the combination = 0.1 is relevant. 

For these values, in the first resonance band we find 


mi 


= 0.995 - A((/)V) - 


The growth index of the first resonance band is thus 


//res ~ (//2 = ^ = 0.0125, 


(2.15) 


(2.16) 


and the following resonance bands are weaker. 

According to (2.15), when M^(t) ~ 1 the first resonance band becomes inactive, and 
the preheating process ceases to be effective. We demonstrate this numerically in Section 3 
in two ways; in the first we insert M^(t) (2.19) by hand; and in the second we simulate Nj 
scalar fields directly as dynamical degrees of freedom. This allows us to test the analytical 
calculation of thermal blocking. 


2.3 Thermal and out-of-equilibrium masses 

The effective mass of the primary particles to leading order in a coupling expansion, is 

-^-^efr(^) = (2.17) 

To this order, the self-interaction from the quartic Higgs coupling is 

MliAt) = A(^V), (2.18) 


with the correlator computed in the out-of-equilibrium state given by the resonance spec¬ 
trum of the preheated Higgs field. Mgeif is therefore time-dependent, and may also have 
a thermal component. The contribution to the effective mass from interactions between 
primary and secondary fields is, to leading order and in equilibrium at temperature T, 

1,2 AT rp 2 

Ml{t) = h\e) = (2.19) 

where the last equality applies for massless helds in thermal equilibrium^. Due to Mgeif and 
Mg, there can therefore be thermal blocking of the resonance in the primary field. This 

^We note that at zero temperature this is probably not a good approximation, since the will have a 
mass proportional to the Higgs vev, At finite temperature, there will also be thermal masses 

for the as well, which at LO and in the high-temperature limit are given by ~ ^ ^ . 
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occurs if either the temperature T or coupling Nfh? are large enough such that the sum 
of mass components (2.17) becomes larger than the driving frequency 

The total effect of interactions with the secondary fields (and self-interactions) is not 
always well described by (2.17) (this is the expression used in [4] for analytical calculations). 
In general, one should instead use the complex, non-local selfenergy S, which appears in 
the corresponding evolution equation in momentum space as 


d} 


+ ikr- 


2 , 2 2 / 
n +g a ' 


[ /s(k,k',t,t')</>k'(tO = 0. (2.20) 

-I Jk' Jt' 


Under certain conditions, Re[S] can be thought of as a thermal mass contribution and 
Im[S] as decay or a damping rate [13]. Such an identification requires that the self-energy 
is approximately local in time. At LO in a loop expansion, S is real, local and equivalent 
to the expression used in the analytical calculations (2.19). However, the complete self¬ 
energy requires more complicated methods to compute, which we implement numerically 
in a classical (as opposed to fully quantum) theory. This classical-statistical approximation 
is known to be very good for high temperatures and large particle numbers [11, 14]. 


2.4 Time scales 

We may identify a number of time-scales in the problem. The basic time-scale is the 
oscillation of the driving field cr, which is Tdriving — "1^^- The resonance builds up in 
a time of the order Tresonance = (10 — 100) Tdriving- The primary particles produced by 
the resonance then decay. If this decay is into other Higgs modes, then the timescale is 
given by the perturbative decay width of the Higgs at finite temperature, giving Tdecay,i oc 

oc If the decay is into excitations of the secondary fields the timescale is 

given by another perturbative decay width, Tdecay ,2 ot T^^ oc . These secondary 

particles then thermalise, and the time-scale for kinetic thermalisation (redistribution of 
the particles into a semi-thermal spectrum) is of order a hundred times the decay time 
''"kin, 1/2 — 100 rdecay,i/ 2 ) ^nd the timescale for complete chemical equilibrium is one or two 
orders of magnitude larger than that, Tchem,i/2 — (10 “ 100) rkin^i/2. In addition, in the 
case of an expanding Universe, there is the Hubble time, th = H~^. 

We can find the ordering of the timescales using a few facts: a succesful resonance 
requires rria > mn', h? = 0(1) (Standard Model couplings); and A = m'jj/(2v‘^) = 0.13. 
We therefore expect that rdecay ,2 ^ ''‘decay,! ^ "^^ 11 , 1/2 ^ ''‘chem,i/ 2 - We also assume that 
the Hubble expansion is negligible, an assumption that depends strongly on what happens 
in the inflaton sector^. The simulations performed here end before kinetic equilibrium has 
completed, so before T]^^in,i/ 2 - 

In the case where nia S> mn and 3> T, the resonance completes much faster than 
any other processes can react. It is then not obvious whether the microscopic processes 
giving rise to the thermal masses can “keep up”. However, the resonance bands are then 

®If the energy is still in the inflaton, the expansion rate could be substantial. If it has been dumped 
into the thermal background of secondary and Higgs fields, they will either be at 100 GeV temperatures, 
in which case expansion is negligible; or they will have been reheated to very large temperatures indeed, in 
which case all mass-scales in the above time-scale estimates must be replaced by that temperature. 
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high above the scales niH and T anyway, and we would expect no thermal blocking, except 
for a small shift in the peak location. We do not discuss this case further. 

However, if ruo- > mn and > T (but not ^), Tresonance — '^decay,i/ 2 ) thermal 
blocking may kick in, if the details of the dynamics generate a mass, and if the excitations 
do not decay away too fast. It is this regime that we investigate in this paper. It is also 
the regime where it is likely that the LO result (2.19) is incomplete. 

2.5 Numerical implementation 

We discretize the system on a 64^-site spatial lattice and solve the classical equations of 
motion in real-time with a time-step of dt = 0.05. All dimensionful quantities will be 
quoted in units of which in lattice units is unity, am^ = 1. This means that the basic 
frequency is well inside the available dynamical range of the lattice, and since we are mostly 
interested in IR physics, the rather large lattice spacing will not be crucial. We generate 
= 64 random classical realisations^ of the 4> and ^ fields for each set of parameters h, T, 
Nf. We run the simulations until t = 1000 which is easily enough to have a strong 
resonance, but shorter than the thermalisation time of the system. Hence we expect to 
see a clear peak in the spectrum, although decays will also have an effect, because the 
time-scales involved are around T^ecay,i /2 

The instantaneous particle number distribution can then be computed from the inver¬ 
sion of (2.8) 

= ((<(>l(k)'/>*(k))(vr|(k)7ri(k))) . (2.21) 


3 Numerical Results 

In the following, we present a number of different approximations to the effective mass. 
We begin by showing the effect of adding a “by hand” thermal mass (section 3.1). We 
then compare this to the case where the secondary fields are dynamical, and show how the 
results depend on the number of fields involved (section 3.2). 

3.1 The LO thermal mass by hand 

Fig. 1 (left) shows the particle number distribution as a function of |kp at time t = 
1000 when the interaction with the fields is turned off, but the self-interaction is 
still present. Note that the x-axis is the momentum squared, and the y-axis is logarithmic. 
We have taken a reference temperature (T/mo-)^ = 2. We see a clear primary peak at 
|k|^/m^ = 0.75 (with some additional structure), on the background of the initial thermal 
spectrum. As expected, it is close but not exactly at the first resonance |k|^/m^ ~ 1, 
because the self-interaction provides an effective mass. There are also secondary peaks 
due to the self-interactions, but it is not possible to see the second resonance band around 
|k|^/m^ ~ 4. The inset figure shows the energy in the Higgs field, which rises monotonically, 

“^In addition to summing over random realisations, we sum over the four Higgs degrees of freedom, and 
over hyper cubic symmetry, i.e. identifying modes of momentum k that only differ by permutations of 
iikx,±ky,±ikz. This gives up to 48 different combinations that are averaged over. 
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Figure 1: The particle number after preheating at t = 1000 corresponding to ap¬ 
proximately 160 inflaton oscillations. Notice the logarithmic scale. Inset is the energy in 
the preheated field(s). The Higgs field is self-interacting and coupled to the “by-hand” 
inflaton, but has no coupling to any other fields. Shown without an additional mass (left), 
and with a mass of = 0.5m^ (right). Note how the additional mass shifts the resonance 
peak to smaller |k|^. 



^ ^ ^ ^ 
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Corresponding M (LO) 


Figure 2: The position of the resonance peaks as a function of the corresponding effective 
(LO) mass, as a function of temperature T (left) and coupling h (right). Left: Black circles 
(top) show only the effect of Higgs self-interaction; red squares (bottom) also includes the 
LO mass; dashed line is explained in the text. Right; The dashed line here is the result 
including only self interaction; points include the LO mass, with fixed = 2 for the Higgs 
field. 


athough not exactly linearly with time. The net transferred energy is not significantly 
affected by the shift in the peak. 

In Fig.l (right) we show the spectrum when we have added by hand a “thermal mass” 
M^(r) = 0.5. The peak is now shifted to = 0.31, which, within the width of the 

peak, is fairly consistent with a shift of 0.5. We still do not have any dynamical secondary 
fields ^ present; their effect is mimicked by the by-hand thermal mass. 

Fig. 2 (left) shows the peak location as a function of the “by-hand” M^(T) of the 
system. To obtain the figure, we repeated the calculation used to make Fig. 1 (left) with 
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varying M‘^{T). The error bars give an estimate of the systematic uncertainty due to the 
width of the resonance peaks in the numerics (compare to Fig. 1). The circular black points 
(top) are equivalent to Fig. 1 (left), i.e. include only self-interactions, calculated at various 
initial temperatures of the Higgs field. As expected, just changing this initial temperature 
changes the location of the peak. The square red points (bottom) are equivalent to Fig.l 
(right) i.e. include a by-hand thermal mass, in addition to the self-interaction. To obtain 
these points, we again varied the initial temperature, but in addition changed M^(T) by 
varying T, while keeping Nfh?, in accordance with (2.19). Thus, the shift represented by 
the square red points is the sum of the effect of the explicit M^(T), and the self-interaction 
of the Higgs, the circular black points. 

The dashed line in the figure shows M^(T) = |k|^. If the effects of thermal masses 
were strictly additive (i.e. as assumed in an analytical calculation), then the position of 
the peak should reach |kp = 0 at the value of M^(T) where the dashed line crosses the 
circular black points, i.e. M^(T) = 0.7. We see that in fact the peak does not leave the 
spectrum until slightly larger M^(T), 0.8 instead of 0.7. Thus, the thermal masses are not 
strictly additive. 

In Fig. 2 (right), we show the location of the peak as a function of M^(T), when we 
imagine that the variation is due to a change in with NfT‘^ kept fixed. In this case, the 
contribution from the self-interaction is constant (symbolised by the dashed line; see Fig. 1), 
and the effect of the explicit M^(T) adds to this contribution. Again, with strictly additive 
behaviour, one would expect that the peak leaves the spectrum at M^(T) = 0.75. There is 
however a small delay to this, and the peak leaves the spectrum at around M^(T) = 0.83. 

In this section, we have shown that a large-enough by-hand thermal mass (2.19) causes 
the resonance to move out of the available modes (mass becomes too big), and preheating 
to stop. The behaviour is convincingly linear, and different effects approximately (but not 
exactly) add up. Higher resonance modes do not seem to provide a sufficiently fast decay 
to prevent this thermal blocking. 

3.2 Dynamical secondary fields 

We now consider the full field dynamics, beyond LO, and investigate whether the resonance 
peak is still shifted out of the spectrum, i.e. we determine to what extent the thermal mass 
is well represented by the LO result. In order to compare to the LO result. Fig. 3 shows 
again the self-interacting (circular black), and LO (filled red squares) results, copied directly 
from Fig. 2 (left). These were the results at LO obtained by varying T. In addition, we 
show two cases with dynamical fields The first (open symbols; second from top) has 
Nf = 1 and a coupling of h? = 3. The second (half-shaded symbols; second from bottom) 
has Nf = 6 and = 0.5. At LO, these results would be equivalent, because N jh?' = 3 in 
both cases. However, beyond LO the degeneracy is lifted. We expect that the large Nf, 
small coupling h? case should be closer to the LO result, since the expression is precisely 
leading order in a coupling and/or 1/A^/ expansion. 

These expectations are confirmed by the simulations. Thermal blocking is still effective, 
but somewhat less so for the fully dynamical system. By the time the resonance is fully 
blocked with the LO mass, the peak in the full dynamics is still at k'^ = 0.4 (Nf = 1) or 
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Figure 3: Shows the position of the resonance peak when varying T, against the corre¬ 
sponding effective (LO) mass (2.19). Filled symbols: Self-interactions only (circles; top) 
and self-interactions plus LO effective mass (squares; bottom). Open symbols (second from 
top): Full dynamical light fields with Nf = 1. Shaded symbols: Full dynamical light fields 
with Nf = 6. 


k'^ = 0.3 {Nf = 6). We also observe that the larger Nf case is indeed closer to the LO 
result. 

Thermal blocking is complete at Nf = 6 for ~ 7.2, whereas for Nf = 1, simple 
extrapolation suggests a of between twice and three times that is required. We dare 
not simulate to larger temperatures, for fear of populating modes near the lattice cutoff. 

We now investigate how the results depend on h, showing the results in Fig. 4. For 
comparison, we show the results including a by-hand thermal mass, keeping = 2 fixed 
but varying h (filled blue diamonds; equivalent to the points in Fig. 2, right). In analogy 
with Fig. 3, we show the corresponding curves including the effects of dynamical fields. 
Open symbols denote Nf = 1 and half-shaded points Nf = 6. In this case, we see that the 
three curves agree well for small couplings h? < 2 {M‘^{T) < 0.3). The dynamical result 
with Nf = 1 bends off at larger couplings, making the LO thermal mass a less robust 
approximation in that case. Even for Nf = 6, the agreement becomes less convincing 
around = 4 {M‘^{T) < 0.6). 

For iVj = 6, it is clear that thermal blocking will complete around M‘^{T) ~ 2, a 
factor of between two and three larger than the LO expectation. For Nf = 1, the picture is 
less clear. However, a linear extrapolation suggests complete blocking around M‘^{T) = 5, 
safely within a factor of ten of the LO estimate. 
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Figure 4: Shows the position of the resonance peak when varying h, as a function of the 
corresponding effective (LO) mass (again through (2.19)). Dashed line: Self-interactions 
only. Filled symbols (bottom): Self-interactions plus effective mass. Open symbols (top): 
Full dynamical light fields with Nj = 1. Shaded symbols (middle): Full dynamical light 
fields with Nf = 6. 


4 Conclusion 

Resonant preheating occurs when the frequency of the driving field (a) corresponds to 
the frequency of some mode of the resonating field ((/>). In many cases, only the first 
resonance mode is effectively resonating. If the effective mass of the resonant field is 
higher than the driving frequency, then there is no preheating because that mode does not 
exist in the spectrum. Thus, the formation of an effective mass is a mechanism that may 
end or prevent preheating. The effective mass can be generated through self-interactions 
or through thermal effects, if the resonant field is coupled to thermalised (or even just 
populated) fields. 

One may expect that this thermal blocking can be modelled by computing the thermal 
mass to leading order in perturbation theory, and that higher order terms contribute only 
small corrections. In this paper we investigated this assumption and found that naively 
using the LO thermal mass is a fair approximation when the number of secondary fields is 
large (> 6). In this case, thermal blocking prevents preheating as the resonating mode is 
shifted out of the spectrum, as represented by 

k^pe.k^k^pe.k(T = 0)-cNfh^T^. (4.1) 

The (model-dependent) constant is c = 1/12 at LO, but is somewhat smaller in the exact 
case (by a factor of 1/2 to 1/10; smaller with increasing Nf). However, when the number 
of secondary fields is small (Nf = 1,2) the LO approximation is poor, and even for large 
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couplings and temperatures, the thermal shifting is ineffective and does not block the 
resonance. 

For models where the primary field is the Higgs, such as the MCH model, the number 
of secondary degrees of freedom coupled to the Higgs field is about 100. Provided that 
it is sufficient to model the fermions and gauge fields by a set of massless scalars, our 
results show that the LO thermal mass gives a good description of the blocking mechanism, 
conservatively within a factor of 2 in couplings and/or temperature. 

Our results are applicable to the thermal blocking of curvaton preheating [4]. In 
this model, the thermal background is generated by the decay of the inflaton, before the 
start of the curvaton resonance. Our results confirm previous analytical results, which 
showed that this thermal background produced substantially blocks curvaton preheating 
[4]. Without this blocking, the MCH model is not viable [5]. There may also be implications 
for preheating in the Higgs Inflation model [15] and for warm inflation [16]. 

A number of unresolved factors relating to the expansion of the Universe have not yet 
been included in our lattice simulations. These include the dilution of decay products, the 
thermalisation mechanism of the particles, the redshifting of modes in and out of bands and 
the decrease of the inflaton oscillation amplitude. In addition, there are outstanding issues 
relating to the field dynamics, including the back-reaction on the inflaton as a dynamical 
field, the decay of primary excitations into secondary fields, and the origin of the initial 
thermal background. All of these await further work, and details will have quantitative (but 
probably not qualitative) impact on the phenomenon of thermal blocking of preheating. 
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